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Abstract: We perform a first-principles calculation of the quantum-limited 
laser linewidth, testing the predictions of recently developed theories of 
the laser linewidth based on fiuctuations about the known steady-state laser 
solutions against traditional forms of the Schawlow-Townes linewidth. The 
numerical study is based on finite-difference time-domain simulations of the 
semiclassical Maxwell-Bloch lasing equations, augmented with Langevin 
force terms, and includes the effects of dispersion, losses due to the open 
boundary of the laser cavity, and non-linear coupling between the amplitude 
and phase fiuctuations (a factor). We find quantitative agreement between 
the numerical results and the predictions of the noisy steady-state ab initio 
laser theory (N-SALT), both in the variation of the linewidth with output 
power, as well as the emergence of side-peaks due to relaxation oscillations. 
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1. Introduction 

The most important property of lasers not captured by semiclassical theories, which treat the 
fields via Maxwell’s equations, is the intrinsic laser linewidth due to quantum fluctuations. 
Above the laser threshold these fluctuations cause a diffusion in the phase of the emitted laser 
signal, leading to a broadening of the observed line, which would have zero width within semi¬ 
classical theory. The magnitude of this linewidth depends upon the geometry of the laser cav¬ 
ity as well as upon the output power of the laser, and was first calculated by Schawlow and 
Townes m, and the standard formula arising from their work, the “Schawlow-Townes” (ST) 
linewidth, is 

= (!) 

where oto is the central frequency of the emitted laser light, jc is the decay rate of the passive 
cavity resonance corresponding to the laser mode, and P is the output power. (Schwalow and 
Townes actually found twice this value in their original work, which assumed the laser was near 
threshold, but it was quickly recognized that far above threshold only the phase fluctuations 
were important, reducing the linewidth by a factor of two). In subsequent decades, improved 
theoretical analyses allowed for the discovery of four significant corrections to this formula. The 
a factor arises from the coupling between intensity and phase fluctuations, and takes different 
forms depending on the nature of the gain medium. For atomic media it was first recognized by 
Lax El and tends to be relatively small; for semiconductor media its importance was realized by 
Henry EHH. and in this context it typically dominates the direct phase fluctuation terms found 
by Schawlow and Townes, and is called the Henry a factor. A second correction arises from 
relaxing the assumption of complete inversion of the gain medium used by Schawlow-Townes; 




this incomplete inversion factor accounts for the actual number of inverted gain atoms El. A 
third and particularly interesting correction is the Petermann factor, which describes the effect 
of the openness of the cavity and the consequent non-orthogonality of the lasing modes M. 
Each of these corrections increases the linewidth from the ST value, which it is natural to 
regard as a lower bound. However, there exists a fourth correction, often referred to as the 
“bad-cavity” factor, which leads to a reduction in the laser linewidth. This correction is only 
appreciable when the cavity decay rate, 70 is fast enough (Q is low enough) that 70 ~ y±, 
where the latter is the dephasing rate of the polarization of the gain medium, which determines 
the gain bandwidth EIEMM). Hence it is appreciable when gain dispersion is significant. 
This correction was first interpreted as a slowing of phase diffusion due to atomic memory 
effects Enani , and subsequently an alternative interpretation was pointed out; the dispersion 
reduces the group velocity of the light within the cavity, leading to an increase in the effective 
cavity Q and a narrowing of the ST linewidth due to the active medium ifT^ . More recently, 
superradiant gain media have been proposed as a way of using the bad-cavity factor to achieve 
ultralow linewidth lasers IBHIll. 

However none of the previous linewidth theories have treated fully the space-dependence 
of the electric fields and the non-linear spatial hole-burning effect in lasers, which greatly af¬ 
fects the stimulated and spontaneous emission rates at different points in the cavity. Recently, a 
steady-state ab initio laser theory (SALT) 12014221 has been developed which treats the spatial 
degrees of freedom essentially exactly, even in the case of multimode lasing. For single-mode 
lasing the validity of the theory only requires that Y± ^ ')^|,where ')^| is the non-radiative relax¬ 
ation rate of the lasing transition; for multimode lasing the condition Y\\ ^ A, where A is the 
free spectral range of the resonator, is also required l22l424l . Subsequently ab initio linewidth 
theories, based on fluctuations in the fields around the SALT solutions, have led to general¬ 
ized linewidth formulas which should be more accurate than the original ST linewidth formula 
with the four previous corrections included as independent multiplicative factors, as is typi¬ 
cally done. The first works of this type used a scattering matrix formulation of the quantum 
fluctuations and input-output theory M25II26I . which captured correctly the generalization of the 
Petermann and bad-cavity factors, but not that of the alpha and incomplete inversion factors. 
The linewidth is expressed in terms of the residue of the lasing pole in the scattering matrix and 
leads to analytic formulas in terms of the SALT solutions. Very recently Pick et al. Il27ll have 
derived a more general analytic formula for the linewidth, by applying a coupled mode noise 
analysis to the SALT solutions. This formula agrees with the results of 12511261 but goes beyond 
them to include correctly more general a and incomplete inversions factors. We will refer to 
this generalized theory, which includes noise effects, as N-SALT (SALT plus noise). We believe 
that the N-SALT linewidth formula quantitatively predicts the laser linewidth (far above thresh¬ 
old) including all corrections in an appropriately generalized form, and in that sense represents 
completely the effects of spontaneous emission on the laser linewidth. We test this hypothesis 
in the current work by direct integration of the laser equations with noise. 

Adding Langevin noise to the steady-state lasing solutions for a gain medium of two-level 
atoms was shown in ll27l to lead to a set of non-linear coupled mode equations for the time- 
dependent fluctuations around the SALT steady-state. Evaluation of the noise-averaged held 
correlation functions from these equations gives the N-SALT laser linewidth in the form: 


Scun-salt = 
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where t/^o(x) is the the semiclassical lasing field inside of the cavity found from SALT, normal¬ 
ized such that / V^Qiix = 1 , and the integral is over the cavity region. e(x) is the total dielectric 
function of the passive cavity plus gain medium, assumed here to be homogeneously broad- 










ened two-level atoms, and N 2 {x) and D{x) are the number of excited atoms and the atomic 
inversion respectively (generalization to multi-level, multi-transition atoms is straightforward 
within SALT and N-SALT, see |I28]|23). a is the generalized a factor IZTl , which can be cal¬ 
culated analytically from knowledge of t/r()(x),e(x) ll27l . This formula is derived under the 
conditions that ScOn-salt ^ This equation reduces to the separable corrections discussed 
above in the appropriate limits 1261271 . but show that in general the incomplete inversion, Peter- 
mann, and bad-cavity linewidth corrections cannot be considered independent from each other 
or of the cavity decay rate. 

Here, we test the predictions of the N-SALT linewidth formula against the Schawlow-Townes 
linewidth formula, including all the relevant corrections by directly integrating the laser equa¬ 
tions using the Finite Difference Time Domain (FDTD) method, including the quantum fluc¬ 
tuations using the method proposed by Drummond and Raymer and employing the time¬ 
stepping method proposed by Bidegaray ED . Many previous numerical studies of spontaneous 
emission in laser cavities have implemented the noise based on knowledge of the lasing mode 
structure E2l - l35l . However, these studies did not have access to the above-threshold lasing¬ 
mode profiles, which can differ significantly from the passive cavity modes used e.g. in cal¬ 
culating the traditional Petermann factor. In our approach we will not make a particular modal 
ansatz. Hofmann and Hess derived FDTD-based noisy lasing equations similar to ours for appli¬ 
cations to semiconductors, but the analysis made further assumptions not valid above the lasing 
threshold ES- The effects of fluctuations in the electromagnetic fields due to thermal noise has 
also been previously studied using the FDTD algorithm llJ7H39l : these effects are necessary 
to include when studying the noise properties of masers or other long wavelength lasers, but 
can be safely neglected at optical frequencies, where the spontaneous emission events being 
considered here dominate the noise of the laser. The approach used in this manuscript is similar 
to that used by Andreasen et al. Il40ti42l . both in the equations used and in the analytic method 
to extract the signal’s linewidth. However unlike those earlier studies Il40l442l we will analyze 
the linewidth far above threshold where it can be compared quantitatively to previous proposed 
formulas. To our knowledge this is the first study of this type. To this end, we will be consid¬ 
ering relatively simple and small laser cavities, allowing us to achieve the spectral resolution 
necessary to resolve the narrow laser linewidths far above the lasing threshold. 

The outline of the remainder of this paper is as follows. In Sec.|2]we demonstrate the equiv¬ 
alence of the macroscopic picture of the N-SALT linewidth formula with the microscopic pic¬ 
ture used by Drummond and Raymer. In Sec.[3]we review the equations and numerical method 
used in the FDTD algorithm to simulate a noisy gain medium coupled to a laser cavity. Sec.|4] 
presents the methodologies for extracting a linewidth from the resultant noisy signal in both 
the frequency and time domains. The results of our study are given in Sec.|5] including the di¬ 
rect comparison between the Schawlow-Townes and N-SALT linewidth predictions in a simple 
laser cavity with single mode lasing, in single mode lasers with a relatively large a factor, and 
in the case of multimode lasing, particularly near the second lasing threshold. Summary and 
concluding remarks are given in Sec.|6] 

2. Microscopic and macroscopic noise equivalence 

There are two different ways of incorporating the effects of spontaneous emission on the electric 
held inside of the laser cavity, either by using the fluctuation-dissipation theorem alongside the 
wave equation, or by including spontaneous emission in the atomic degrees of freedom, which 
are coupled non-linearly to the wave equation. In this section we will explicitly demonstrate the 
equivalence of these two methods, which we term the macroscopic and microscopic perspec¬ 
tives respectively, as the derivation of the N-SALT linewidth equation uses the former method, 
while the Langevin equations augmenting the FDTD simulations use the latter. This section 




also serves as a proof that despite the non-equilibrium nature of the laser, with power flowing 
in and light flowing out, the system does reach a point of stability wherein the fluctuations of 
the electric held can be appropriately treated with the fluctuation-dissipation theorem. 

The derivation of the N-SALT equation incorporates all of the noise due to the quantum 
fluctuations in the gain medium directly into the wave equation as m 

[V X V X -cu^e(a),Eo)] E = (e(a),E) - e(a),Eo))E-f Fj, (3) 

where e(a),E) is the full dielectric function of the cavity and gain medium, e(a),Eo) is the 
non-linear saturated dielectric function of the cavity evaluated using the semiclassical lasing 
mode Eo(x) = v/7v^o(x), where I is the lasing mode intensity, and F 5 is a random noise source 
corresponding to the spontaneous emission from the gain medium. The first term on the right 
hand side of Eq. Q corresponds to the effective source due to fluctuations in the held leading 
to fluctuations in the saturation of the gain medium, while the second term corresponds to 
spontaneous emission contributing directly to noise in the electric held. The inclusion of the 
full space-dependent non-linearity of the active cavity dielectric function above threshold in 
the noise term is a key feature distinguishing N-SALT from previous linewidth theories. The 
autocorrelation of the random noise source is assumed to be given directly by the fluctuation- 
dissipation theorem, 

(Fj(x,a))Fs(x',a)')) = 2/ttu"^Im[e(tu,Eo)]coth 5{x-x')5{(0- (o'), (4) 


where i3(x) = (l//!aJo)ln(Ai(x)/A 2 (x)) is the effective (negative) inverse temperature of the 
inverted gain medium, with Ni and N 2 are the number of atoms in the ground and excited atomic 
levels respectively. (Note that Im[e(a),Eo)] < 0 in the inverted state, so that the correlation 
remains positive). 

In this treatment of the noise in the laser held due to spontaneous emission, the atomic de¬ 
grees of freedom have been completely integrated out, and the fluctuation-dissipation theorem 
has been invoked from a macroscopic perspective, relating the autocorrelation of the noise 
source to the imaginary part of the material response function and a temperature dependent 
term. The hyperbolic cotangent factor arrises as a sum of a Bose-Einstein distribution and a 
factor of 1 /2 from the quantum zero-point fluctuations, which is why the auto-correlation does 
not vanish in the zero temperature limit (j3 —?► °°). However, it was shown by Henry and Kazari- 
nov that the contributions from the zero-point fluctuations cancel in the linewidth formula Il43l 
(a simpler, semiclassical proof of this is in Ref. IZTl ). and as such it is convenient to explic¬ 
itly subtract this contribution, allowing for the effective temperature of the gain medium to be 
determined by relative occupations of the atomic levels comprising the lasing transition. 
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where D{x) = A 2 (x) — Ai(x) is the number of inverted atoms. Thus, for the laser systems 
considered here, Eq. (0 can be written as 


(F^(x, a))F 5 (x', to')) = 4 r 2 a)"^Im[e(x,«)] 



5{x — x')5{(0 — (o'). ( 6 ) 


In contrast to this macroscopic picture, many traditional theories of the noise due to spon¬ 
taneous emission from the gain media begin by treating the Langevin forces on the quantum 
operators of individual gain atoms and building up an understanding of the total noise this gen¬ 
erates in the electric held, a more microscopic viewpoint ||2l[30ll44l. We will demonstrate the 










equivalence of these two methods by deriving the total Langevin force on the polarization from 
the microscopic perspective. For a two-level atomic gain medium, the evolution equation for 
the off-diagonal matrix element of the ath atom, p^j , including the Langevin force, Fj^j (f), is 
given by, 

a,p(“) (f) = - (7x + ^ ® , t)+ ), (7) 

in which a)„ is the atomic transition frequency, 7 x is the dephasing rate, and 6 is the dipole cou¬ 
pling matrix element. Furthermore, the evolution of the inversion for that atom, including 
the Langevin force, (f), is given by 


- £/(“)) + Ae. E(x(‘= 


,f)(p(“)*-p(“')+r: 
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( 8 ) 


where is the inversion of the ath atom in the absence of any electric field. Finally, the wave 
equation for the electric field can be written in this context by explicitly including the coupling 
between the field and each individual gain atom (see Eqs. (5.48) and (5.55) in Ref. 1041), 

[V X V X -(oIec] E(x,a)) = 47 ra)o 6 ^ 5 (x-x(“*)p 2 “*, (9) 
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in which we have approximated that the electric field is oscillating at frequencies close to the 
semiclassical lasing frequency, coo, and retained only the positive frequency components for 
both the electric field and atomic polarization. Our aim is to determine the form of the effective 
total Langevin force on the electric field by solving Eqs. (|7]i and ® for the polarization and 
inversion, insert these expressions into the wave equation, and collect the resulting Langevin 
force terms. 

To leading order, p 2 i will oscillate at the lasing frequency, (Oq, and if we approximate this as 
its only frequency component, we can solve for 


{“) 

P21 = 




h{(Oo-(Oa + iY±) 
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where the electric field is assumed to be a constant over the volume of the atom at x(“). The 
fluctuation dissipation theorem states that the strength of the fluctuations is proportional to the 
strength of the dissipative terms. Thus, for the Class A and B lasers considered here, Y\\ < 7±> 

so <C r|pj(f), and we can safely ignore the fluctuations in the atomic inversion. Thus, 

we can insert Eq. (fTOl i into Eq. (|9]l. 
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Equation (fTTT i allows for the identification of the spontaneous noise in the polarization, Pat, 
using Eq. Q and noting that F 5 = —AkgP'En, as 


PAr(x,tU) = ^5(x-x(“)) 




(OQ-(Oa + iYl_ 




( 12 ) 










We can now directly calculate the correlation function of the spontaneous noise in the polariza¬ 
tion using the correlation of the atomic Langevin force M, 




in which each atom is taken to only be in equilibrium with its reservoir Eol . Note that the only 
place the non-equilibrium nature of the reservoir comes in is via the term — p“ . 

Since no higher moments or correlations enter the calculation, it is safe to dehne an effective 
temperature for the system which can be negative via this relation and apply the fluctuation- 


dissipation theorem. Note also that contains the non-linear effect of gain saturation and 


spectral hole burning when calculated by the FDTD method given below. 

By assuming that the inversion is relatively stationary, we can identify the same frequency 
auto-correlation of the noise as na 




(14) 


in which we have again dropped the noise source proportional to to be consistent with the 
approximation neglecting fluctuations in the inversion made above. This allows us to solve for 



(15) 


where the number of atoms in the upper lasing state, N 2 (x) has been identified using. 


^ rv 


(16) 


a 


Upon substitution of the imaginary part of the dielectric function. 



(17) 


we can identify the same frequency auto-correlation of the noise source F 5 as 



(18) 


Finally, noting that the different frequency auto-correlation function can be found as m. 



(19) 


and using the definition of the temperature factor given in Eq. (|5]l, we recover the expected 
auto-correlation of the random noise source given in Eq. (| 6 ]l. With this, we have verified that the 
microscopic and macroscopic methods of treating the fluctuations in the gain medium produce 
identical results, which allows us to use a microscopic model of the gain medium in our EDTD 
simulations to test the predictions of the N-SALT theory. 


3. FDTD equations 

Having now demonstrated the equivalence of the microscopic and macroscopic fluctuation 
models, in this section we show how to include the microscopic fluctuations of the gain 







medium within an FDTD simulation of a laser. The FDTD algorithm has been known since 
the 1960s ll46l and is ubiquitous across many fields of study El. However, only a few pre¬ 
vious works have used the algorithm to study the noise in lasers ll40H42l . and none (to our 
knowledge) have studied the linewidth far above the lasing threshold as we do here. For this 
reason we will briefly review the simulated equations here. The Maxwell-Bloch equations for 
a two level atomic gain medium in a one dimensional cavity can be written as 
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where E„ and B„ are the electric and magnetic field densities at the spatial location Xn within 
the lasing cavity, Vq is the volume associated with each grid point, 7^ is the total atomic off- 
diagonal density matrix element (related to the polarization) with a positive frequency compo¬ 
nent, D„ is the inversion of the N„ atoms at the spatial location Xn, is the inversion in the 
absence of an electric field and plays the role of the effective pump strength in this theory, and 
py'^ and Fn^'^ are the Langevin forces experienced by the atomic off-diagonal density matrix 
element and inversion respectively. The choice of 7„ for the total off-diagonal density matrix 
element is made for ease of comparison with Drummond and Raymer, who use 7,7 to denote 
the same quantity, and is defined as 


W = P 2 I ^ = ^"P21 W ■ 


(24) 


The Langevin forces can be written as 
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in which 721 is the pumping rate from lower level 11 ) to | 2 ) and is given by. 


721,« = 


n 

2 


Do^n 

IP 


(27) 


and 7 /> = 7 x — /2 is the pure dephasing rate. Randomness is introduced through the stochastic 

variables ^, which are complex except for G K, and satisfy ll^ 


{p{t)yj\t')) = 5{t-t')5„n,5ij. (28) 

Many of the terms in Eqs. (I25l426l l stem from resolving the dilemma of the operator order¬ 
ing when reducing operator equations to c-number equations. However, for studying the laser 
linewidth above threshold, the difference caused by this ambiguity is minimal, as the addition 














or removal of a spontaneous emission event is negligible in the presence of the large number of 
gain atoms necessary for lasing to occur. Thus most of these terms are expected to be negligible, 
an assumption which we check a posteriori after retaining the leading terms: 


py^ =yy ( 29 ) 

py ( 30 ) 

Finally, in accordance with the discussion in the previous section, the thermal fluctuations of 
the electric and magnetic fields have been neglected. 

The Maxwell-Bloch equations can then be discretized for use in the FDTD algorithm fol¬ 
lowing the weak coupling method proposed by Bidegaray lISTl . evolving the atomic variables 
simultaneously with the magnetic field, but at the same spatial locations as the electric field 
so as to avoid solving a non-linear equation. Furthermore, it is useful to separate the real and 
imaginary components of the atomic off-diagonal density matrix element, = y -i-ij^, 
resulting in 



where Un = is the vector of the atomic variables, d,, = {Y\\Do.n,0,0) is the pump¬ 

ing vector, 1 is the 3x3 identity matrix, M is a matrix which contains the coupling information 
between the atomic variables, 

/ _^| 0 -fPniti) \ 

M=\ 0 -7x (Oa , (34) 

\ —(Oa —Y± / 

and f„ is the Langevin force vector, whose elements are 



/«,2 ^yp{P>n{t,_l^)+N„) + y'> yyynNn, (36) 


fn3 ^yp{P>n{t,_l)+Nn) + s/yynNn- (37) 

where we have renumbered the random variables y'\ which continue to satisfy Eq. (l28T l. but 
are now real, rather than complex, and introduced a factor of in this conversion process 
(except for y \ which was real to begin with). Here we have used the final approximation that 
the Langevin force vector only depends upon the inversion at the previous time step, rather 
than the average of the previous and current time steps which would result in a non-linear equa¬ 
tion II 42 I . This is justified for the simulations performed here because the inversion, d„, is many 
















orders of magnitude smaller than the total number of atoms, Nn, and thus these inversion depen¬ 
dent terms will have minimal impact upon the overall strength of the noise. For the discretized 
Langevin forces, the stochastic variables are chosen from a standard uniform distribution, 
and then renormalized to satisfy 

= (38) 

Eas. lSTllJTl can now be readily evaluated numerically. 

4. Linewidth analysis 

Broadly speaking there are two main ways of extracting a linewidth from a noisy signal; by 
either fitting a curve to the frequency domain data or calculating the cross-correlation of the 
time domain data ll48l . Here we will use both methods an compare them; first we calculate 
a linewidth from the spectral data and then confirm this linewidth by calculating (0(f')0(f)), 
where 0(f) is the phase of the electric field. 

4.1. Frequency-domain analysis 

To analyze the spectrum of the electric field output from the cavity, E{co), and find a linewidth, 
we will use the method proposed by Andreasen et al. Il42l . and fit the spectrum to a Lorentzian 
through the use of an error function. We assume that the noise is Lorentzian, 

L{(0)=(—)- -^^-2-2 (39) 

\n J {(0 - (Oq)^ s^ 

where s is the half-width half-maximum of the noise, ScUfdtd = 2s. The Lorentz error function 
can then be defined as 



As such, this integration can be carried out numerically directly upon E{(o), and then fit to 
Eq. (l40t . Eor all of the data shown in this paper the curve fitting is carried out using itera¬ 
tive least squares estimation. Performing this integration requires knowledge of the lasing fre¬ 
quency, CUo, which is known from the semiclassical SALT calculation. However, the presence 
of noise results in a slight shift of the semiclassical lasing frequency IJI, and the slightly differ¬ 
ent discretization schemes used between the SALT and LDTD calculations yield an additional 
shift in the lasing frequency, which together lead to a slightly shifted integrated spectrum, both 
horizontally and vertically. As such it is useful to include two other fitting parameters in the 
Lorentz error function. 
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where d plays the role of the horizontal offset and c is the vertical offset. Using this correction, 
the calculated linewidths are robust to the choice of dp so long as the curve fitting algorithm 
converges. 

An example of this process can be seen in Lig. [T] where the left panel shows the spectrum 
of the output electric field for a dielectric slab cavity. To compute the power spectrum, or tech¬ 
nically the periodogram 1481 of the noisy signal, we chop the simulated time-domain field 
E{t) into ~ 10 pieces and perform a discrete-time Lourier transform (DTLT) l49l on each 
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Fig. 1. (a) Intensity spectrum of the output electric field of an « = 3 dielectric slab cavity, 
shown in the schematic. The simulation parameters for the cavity are y± = .5, COa = 42.4, 
= .01, 6 = 2 X 10^®, Na = 10*°, and the cavity is uniformly pumped at Dq = 0.275 
which is close to 5 times the threshold lasing pump of Dq iIj^ = 0.0488. The rates quoted 
here are given in units of c/L, while the intensity is given in SALT units of 46^/(ft^yx]^|), 
and the number and inversion of gain atoms are given in the SALT units of 4n6'^/{tjY±). 
(b) Plot of the fitted Lorentz error function (red line) and numerically integrated FDTD 
data (blue dots) of the simulation shown in (a). The spectral resolution for the simulated 
data in (a) and (b) is dco = 1.96 x 10^^. The analytic curve fit parameters are found using 
MATLAB’s curve fitting algorithms. 


constituent piece, and then ensemble-average the resulting spectra \E{co)\^ using Bartlett’s 
method HSll . The right panel shows the Lorentz error function integral calculated numeri¬ 
cally and fit against the analytic curve. The resulting linewidth predicted by this method is 
5(Bfdtd = 2.22 X 10^^, which is around an order of magnitude larger than the resolution of the 
resultant spectra, d(0= 1.96 x 10^^, given in units of c/L. 

4.2. Time-domain confirmation 

This calculation can be independently confirmed by calculating the autocorrelation of the output 
electric field as a function of time and expressing this as a function of the phase correlation, 
which is defined in terms of the linewidth of the signal. Writing the output electric field as 

£(f) =Ccos(a)f-F0(f)), (42) 

The autocorrelation of the electric field, REsiSt) = {E{t -\- 5t)E (f)), can then be written as 

REsiSt) = {E^{t)cos{co5t5(j){5t))) — {—sm{2cot-\-2(j){t))sm{co5t-\-5<j){5t))), (43) 

where the double angle formula has been used in finding the second term on the right hand side, 
and 50 (5f) = 0(f -F 5t) — 0(f). By assuming that the phase shift 50 (5f) is uncorrelated with 
the phase 0(f), we can separate the correlations, note that the second term averages to zero, and 
again apply a trigonometric identity, resulting in 

^££(5f) = — [cos(tu5f)(cos(50(5f))) — sin(a)5f)(sin(50(5f)))]. (44) 

This assumption that the phase shift, 50 is uncorrelated with the instantaneous phase, 0, is 
analogous to assuming that the gain medium has no memory effect, and is consistent with the 
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Fig. 2. Plot of the autocorrelation of the electric field simulated numerically for the same 
parameter used in Fig. [T] (hlue line) and the analytic prediction for the envelope of the 
autocorrelation given in the second factor in Eq. l l46b (green line). The fast oscillations in 
the numerically simulated electric field are at the lasing frequency coq, which is much faster 
than the other time scales in the problem and leads to the densely packed curve shown in 
blue. Quantities are normalized, and plotted in units of dcodl. 


earlier assumption that the bad-cavity factor is unity for the systems studied here. The second 
term in Eq. (l44li averages to zero as well, as the phase shift is equally likely to be positive or 
negative. Finally, the cosine of the phase shift can be Taylor expanded, and noting the definition 
of the linewidth, 

{5(j)^{5t)) = 5co5t, (45) 

the electric field autocorrelation can be written as 


REE{5t) = — cos(a)5f) 


5(o5t 


0{5t^) 


(46) 


showing that in the presence of phase diffusion, the correlation should decrease linearly for 
small 5t. 

This trend can be observed in Fig.l^for the same simulation as shown in Fig.[T] where the 
prediction for REE{St) is evaluated using 5(0 found by the frequency domain method from 
the previous section and Eq. (|46] | (green line), and numerically calculated (blue line). The fast 
oscillations seen in the numerical data are due to (o5t ^ 1, and are predicted by the theory 
derived above. The semi-quantitative agreement seen between the frequency domain linewidth 
prediction and the time domain prediction calculated here provides a consistency check, though 
we will use the frequency domain method for the remainder of the calculations performed here. 


5. Results 

To test the predictions of the N-SAFT linewidth, given Eq. (|2|i, with the Schawlow-Townes 
linewidth HI, we first study the simple one-dimensional, single-sided dielectric slab cavity, 
n = 3, used in the previous two sections in Figs.[T]and|2] Here, we use the “fully-corrected” form 
of the Schawlow-Townes linewidth as the point of comparison, which includes the Petermann 
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Fig. 3. (Left panel) Plot showing the linewidth predictions given hy the N-SALT given in 
Eq.dD (green), corrected Schawlow-Townes theory given in Eq. (|47j (blue), integral form 
of the Chong-Stone linewidth formula given in Eq. ( 148b (orange), and EDTD simulations 
(magenta) for a uniformly pumped, dielectric slab cavity with n = 3, COa = 42.4, = 

.5, = . 01 , 6 = 2 X 10 ^^, and = 10 '°. All of the linewidth formulas are evaluated 

using the spatially dependent integral definition of the power given by Eq. (gg. (Right 
panel) Plot of the same data shown on a log-log scale, with reference lines for strict inverse 
power dependence, provided for comparison (black dashed). Schematic inset shows 
the cavity geometry. The rates and frequency are given in units of c/L, the number of atoms 
in the cavity is given in terms of the SALT units of 4;r6^/(ftyx), and the output power is 
given in the SALT units of 48~ 


factor, bad-cavity correction, and Henry a factor, and is given by, 

2 

+ (47) 

where ^o(x) is the passive cavity resonance corresponding to the lasing mode, the spatial av¬ 
erage of the inversion and occupation of the upper lasing state is denoted as D = f D(x)dx, 
the spatially averaged inversion is used to calculate the bad-cavity factor, and a is the Henry 
a factor. The first term in parentheses of Eq. (l47l i corresponds to the cavity-averaged incom¬ 
plete inversion factor and the second corresponds to the Petermann factor ||5]|2^ . The quantities 
t/fo(x), 0o(*)’ £(*) ^re calculated using SALT, while the FDTD linewidths are ex¬ 

tracted using the method described in Sec. 14. II and run for enough time steps to average together 
at least six resulting spectra using Bartlett’s method. For the chosen parameters Yc ~ Y\\ placing 
it on the border between Class A and Class B lasers lISOll . close enough to the former that no 
relaxation oscillation side-peaks are seen in the resulting spectra. 

As can be seen in the left panel of Fig. [2 excellent quantitative agreement is seen between 
the N-SALT prediction (green line) and the linewidths measured through direct integration 
of the noisy Maxwell-Bloch equations (magenta triangles), while both results differ from the 
corrected Schawlow-Townes theory (blue line). This discrepancy is shown to be more than a 
simple scaling factor in the right panel of Fig. [2 where the same data is plotted on a log-log 
scale, and it can be seen that the power law narrowing of the linewidth with respect to the output 
power differs between the N-SALT and corrected Schawlow-Townes linewidth predictions. 
Somewhat surprisingly only the N-SALT and FDTD results are very close to P * (black dashed 
lines), the others are show a measurably faster narrowing. 

To understand the source of this discrepancy, we also plot the Chong-Stone linewidth 


T^coo)^ 

(N2\ 

J\(j>o{x)\^dx 

2 

1 

2P 

[dJ 

J(l>l{x)dx 


1 1 ^ (9e 1 
^ ' le 1 ^ 

















calculated using its integral form ll26l . 



( 48 ) 


where we have neglected the vanishingly small boundary term (see 12^). The Chong-Stone 
linewidth formula is derived through considering the behavior of the SALT-based scattering 
matrix of the cavity, and thus is able to account correctly for all effects stemming from the 
cavity; it gives the proper cavity decay rate above threshold, and the same Petermann factor, 
and bad-cavity correction as N-SALT. However, it does not provide an accurate treatment of 
the fluctuations inside the gain medium, particularly amplitude fluctuations, and thus is unable 
to And the a factor and finds an inaccurate, cavity-averaged incomplete inversion factor similar 
to conventional theories. For the dielectric slab cavity studied here, the detuning of the lasing 
mode from the atomic transition is very small, such that a ^ 1. Thus the significant discrep¬ 
ancy between the N-SALT and FDTD results and the Chong-Stone prediction indicates that 
the largest source of discrepancy lies in the treatment of the incomplete inversion factor. The 
ratio of the N-SALT and Chong-Stone linewidth predictions in the limit that a = a = 0 can be 
written as 


gcucs ^ f-./^(x)|yo(x)P^x 
Scun-salt /A2(x)|i/ro(x)|2£/x 


(49) 


However, for the two-level atomic gain media simulated here, the number of atoms in the ex¬ 
cited atomic level is nearly constant N 2 ~N\~ N j2, allowing for this ratio to be expressed 
as 


5(0cs _ /L>(x)|yo(x)|2^x 


(50) 


5®n-salt /|V'o(x)P'^x/D(x)(fx' 


In absolute terms, the fluctuations in N 2 , N\, and D are all of the same magnitude, but as 
D{x) = N 2 {x)—Ni (x) <C N 2 , its spatial variation is much larger on a relative scale and cannot be 
neglected, leading to a significant discrepancy between the N-SALT/FDTD and Chong-Stone 
linewidth predictions. Note that the approximation of spatial invariance of the occupation of 
the upper lasing level does not necessarily hold when considering more realistic gain media, 
with more than two levels, and is a result of the well known difficulty in pumping a two-level 
medium past the transparency point to achieve lasing. However the residual discrepancy be¬ 
tween Chong-Stone and the corrected ST prediction indicates that the incomplete inversion 
factor only accounts for roughly half the discrepancy, and the remainder (Petermann and bad- 
cavity effects) would be present in lasers with more than two levels. 

The implications of the relation expressed in Eq. (fSOl) can be understood graphically from 
Fig.m where the left panel shows the steady-state inversion, Z)(x), within the cavity for dif¬ 
ferent values of the output power generated by the cavity, and the right panel shows the spatial 
dependence of the lasing mode profile, |v^q(x)|, for the same values of the output power. As 
the pump on the gain medium. Do, is increased, the amplitude of the held within the cavity in¬ 
creases, as does the output power. However, due to spatial hole-burning in the gain medium, the 
impact of the higher held intensity within the cavity is not felt uniformly in the inversion; thus 
the average inversion within the cavity still increases as the pump is ramped, mostly due to the 
positions near the mirror in the cavity where the electric held is very weak, while the weighted 
average of the inversion with the held intensity remains relatively constant, as the inversion 
where the field intensity is maximized stays relatively constant as the pump is increased. Thus 
as noted, we do expect to see the corrected Schawlow-Townes and Chong-Stone linewidth pre¬ 
dictions decrease faster than 1 /P, as is observed in the right panel of Fig. [3l as both the output 
power, P, and spatially averaged inversion, D, increase as the pump strength. Do, is increased 
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Fig. 4. (Left panel) Plot of the steady-state inversion, D{x), as a function of the location 
in the cavity for three different values of the output power, P = 0.524 (blue), P = 1.252 
(green), and P = 3.116 (red). These values correspond to the first, sixth, and eighteenth 
data points shown in Fig. [3 Strong spatial hole-burning is seen in the inversion due to the 
lasing mode. Schematic depicts the cavity from Fig. [3] (Right panel) Plot of the normalized 
spatial prohle of the lasing mode, |t//o(x)|, as a function of position in the cavity for the 
same three values of the output power shown in the left panel. The output power is given 
in dimensionless SALT units of 46^/(/!^7 j^')^|). 


(see Eq. (l48Tl). In contrast, the integral of the inversion weighted against the field intensity, used 
in the N-SALT linewidth prediction, does not change as the pump is increased. Thus, even 
for the two-level atomic gain medium studied here, N-SALT gives a good 1/P line narrowing. 
Siegman has previously suggested that the incomplete inversion factor might lead to deviations 
from the strict inverse dependence of the laser linewidth upon the output power, but was unable 
to test this hypothesis ED. 

We note that it is important in these comparisons to calculate the output power from its 
fundamental definition via Poynting’s theorem E2ll . 

T’= |^/lm[-e(x)]|Eo(x)pr/x, (51) 

where this equation is given in Gaussian units, Eo(x) = is the unnormalized lasing 

mode, and I is the mode intensity. Performing this calculation relies on finding the correct 
space-dependent quantities, which can be obtained using SALT. The quantitative agreement 
seen between the N-SALT linewidth prediction and the LDTD simulations shown in Lig.|Dpro- 
vides independent confirmation that this is the correct formulation of the output power to use. 
However, in many treatments of laser emisson, which do not treat the full space dependence, 
the output power is calculated using Eol 

^ST = Ycnhcoo, (52) 

where n is the average number of photons in the cavity. Using this form of power calculation can 
introduce a substantial error; using the corrected Schawlow-Townes theory with this spatially 
averaged power for the parameters of Lig. [S] leads to a linewidth roughly a factor of two larger 
than the N-SALT and LDTD results. Thus we see that it is critical to use all of the spatial 
information in the fields Eo(x) and D{x) obtained from SALT in order to quantitatively predict 
the laser linewidth. 
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Fig. 5. (Left panel) Plot showing the linewidth predictions given by the N-SALT (green 
line), corrected Schawlow-Townes theory (blue line), and FDTD simulations (red diamonds 
and magenta triangles) for a uniformly pumped, dielectric slab cavity with n = 3, 0 ),, = 42.4, 
7x = -5, = .04, 6 = 4 X 10^®, and = lO**', as shown in the schematic. The results 

of the new FDTD simulations are shown as red diamonds, and are plotted alongside the 
FDTD results from Fig. shown as magenta triangles. (Right panel) Plot showing the 
linewidth predictions given by the N-SALT (green line), rescaled N-SALT prediction from 
Fig. [U (magenta dashed line), corrected Schawlow-Townes theory (blue line), and FDTD 
simulations (cyan squares) for a uniformly pumped, dielectric slab cavity with n = 3, (Oa = 
42.4, yx = .25, )(| = .02, 6 = 2 x 10^®, and Na = 10*°. The rates and frequency are given 
in units of c/L, the number of atoms in the cavity is given in terms of the SALT units of 
4nd^/{hY±), and the output power is given in the SALT units of 46^/{TP'YI-Yw)- 


5.1. linewidth scaling relations 

The overall intensity of the electric field enters directly into the linewidth formulas only through 
the output power, Eq. ( fSTT ). SALT demonstrates that the electric field can be written in terms of 
dimensionless units, Eo(x) = y\\l29 )Esalt(x) II221I53L and thus the output power can 

also be written as, 

P = II /Im[-e(x)]|EsALT(x)|2«fx. (53) 

This is how the dimension-full parameters stemming from the properties of the gain medium 
directly enter into all of the linewidth formulas discussed here. In particular we can rewrite the 
N-SALT linewidth in SALT units as, 

( 402 ^ hcoQ cdo2/Im[e]|vfoP^^x/Im[e]§|i/r(,|2£/x 

oon-salt = 72 - -;- 7 -- 7 —-^t+ « ). ( 54 ) 

2Lsalt |jvr2(e-Ff^|^)«fxp 

where Lsalt is the output power calculated using the electric field measured in SALT units. 

Using SALT units and the stationary inversion approximation implies powerful scaling re¬ 
lations between lasing solutions at different gain medium parameter values M201I22II . Similarly, 
Eq. (|54] | implies various scaling relations for the linewidth. It separates out the dependence of 
the intrinsic laser linewidth upon 0, Y\\, and the leading dependence upon Yj_ and thus predicts 
that the linewidth should obey a set of scaling relations. Lor example maintaining the ratio of 
Y\\/0^ should yield the same linewidth, and keeping the ratio Yi_Y\\/9^ constant should result 
in only very modest changes in the linewidth (changing Yj_ only changes the strength of the 
bad-cavity correction). These predictions are confirmed by LDTD simulations. In the left panel 












of Fig.|5] the linewidth is calculated via FDTD for a different value of and 9 than in Fig. 3, 
while keeping the ratio Y\\/(^^ to that in Fig.|3] The resulting FDTD linewidth (red diamonds, 
plotted alongside magenta triangles from Fig. |3ll is seen to be identical. This serves as a vali¬ 
dation of the FDTD simulations shown here, as 9,Y\\ enter into the equations in a non-trivial 
manner from which the scaling relations are not apparent. 

In practice however checking these scaling relation may be difficult, as the total relaxation 
rate of the inversion, Y\\, can be written as a sum of contributions from spontaneous emission 
and non-radiative decay, 

Y\\ ~ Jspon + Ynr-, (55) 

in which the spontaneous decay rate can be written as 


Yspon — 


AafsCO^nG^ 
3^2 ’ 


(56) 


where a/j is the fine structure constant and Yspon is seen to be exactly dependent upon 9^. 
Thus, in the limit of an atomic gain media without a non-radiative decay channel available 
from the upper level to the ground state, the ratio of Y\\/(^'^ in the linewidth does not yield 
any new information as these two parameters are not independent. However, this analysis does 
verify the intuitive statement that the laser linewidth will be reduced if the non-radiative decay 
rate is substantially larger than the spontaneous emission decay rate, decreasing the overall 
significance of spontaneous emission to the system, as the relative ratio of 9^/Y\\ that appears 
in Eq. (fS^ will be reduced. 

As noted, the scaling of the linewidth with the ratio is not exact as there is an 

additional dependence on Y± in the bad cavity factor once 7c 7±. In the right panel of Fig.|5] 
the ratio of Yi-Y\\l^'^ is held constant and equal to that in the left panel of the figure, but Y± is 
decreased so as to make the bad-cavity factor significantly different from unity. Thus instead of 
remaining constant the linewidth decreases, in this case by roughly a factor of 2/3. However it 
is possible to account for this failure of scaling by including a further approximate scaling by 
noting that when Wa ~ Oo ^ 7^, we can express the bad-cavity factor as 
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Using this, we can rescale the N-SALT linewidth prediction by calculated using the 

simple form on the right-hand side of Eq. (fSTl i (magenta dashed line). With the additional rescal¬ 
ing the N-SALT linewidth for the parameters of the left panel now agrees with the N-SALT 
prediction for the new gain media parameters in the right panel (green line) and quantitatively 
agree with the EDTD simulations (cyan squares). This also verifies that the N-SALT form of 
the bad-cavity factor correctly reduces to previously known approximations ll2l fT^fT3]l . at least 
for the parameters chosen here. 

5.2. relaxation oscillation sidebands 

In Class B lasers, fluctuations in the amplitude of the electric field undergo relaxation oscilla¬ 
tions while decaying to the steady-state. These relaxation oscillations give rise to side-peaks in 
the spectrum of the output intensity and in this section we will demonstrate that the N-SALT is 
able to correctly reproduce the location and size of these side-peaks EH. It has been known for 
many decades that the relaxation oscillation frequency increases as the laser is pumped further 
above threshold ll55l . but previous studies did not take into account the spatial variation in the 
gain saturation, which was shown lIZTll to play an important role in quantitatively predicting 







the laser linewidth in Sec. |5] Using the spatial lasing mode profiles and inversion calculated 
using SALT, N-SALT demonstrates that the output intensity spectrum is dependent upon the 
total local decay rate im, 


which contains contributions from both the non-radiative decay rate of the inversion, , as 
well as the local rate of stimulated emission given by the second term in Eq. (fSSl) . N-SALT 
yields two main results for the effects of relaxation oscillations on the linewidth. First, that 
relaxation oscillation side peaks will appear for cavities whose parameters satisfy the inequality 
5a)N-SALT 'Ki < / A{x)dx, in which 


A(x) = 2/Re 
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where I is the intensity of the electric field, as defined above. 

Second, N-SALT gives an explicit form for the output intensity spectrum in the presence of 
relaxation oscillations (with a = 0): 
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The second term in Eg. [^describes the side peaks due to relaxation oscillations. In Fig.|6] 
we show the output intensity spectrum of a dielectric slab cavity pumped above the first lasing 
threshold, in the parameter regime where side peaks are expected. Each of the plots shows a 
comparison between the N-SALT prediction (red line) and the FDTD simulations (blue line) 
for increasing values of the pump, (a) to (c). As can be seen in all three plots, excellent quan¬ 
titative agreement is seen between the simulated spectrum and the N-SALT prediction. To re¬ 
iterate, N-SALT has no free parameters, so the agreement seen here is a demonstration of a 
first principles test of N-SALT. As can be see in the FDTD simulations, there are additional 
peaks in the spectrum at a distance of twice the relaxation oscillation frequency from the cen¬ 
tral peak. In principle N-SALT can be used to predict these additional side-peaks as well. Fi¬ 
nally, relaxation oscillations are proportional to the square root of the decay rate of the cavity, 
Oro ^ a/(1 /L) f Y(x)dx, thus we expect for the side peaks seen in the spectrum to move away 
from the central peak as the rate of stimulated emission increases due to an increasing pump. 
As the pump is increased from Fig.|6|a) to Fig.|6jc) we observe exactly this behavior in both 
the FDTD simulations and N-SALT results, verifying this prediction. 


5.3. large alpha factor 

The a factor accounts for the phase fluctuations due to changes in the susceptibility of the 
gain medium from intensity fluctuations, and is known to be quite large in semiconductor gain 
material, where it is referred to as the Henry a factor. The N-SALT linewidth theory is quite 
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Fig. 6. Plots showing a comparison between the N-SALT prediction (red) and FDTD simu¬ 
lations (blue) of the intensity spectrum for increasing values of the pump, Dq, for a single¬ 
sided, dielectric slab cavity with n = 1.5, = 40.7, yx = 1, = 0.0025, 6 = 6x 10^*®, 

and Na = 10*®. (a) Dq = 0.18, (b) Dq = 0.28, (c) Dq = 0.38. As can be seen, increasing 
the pump value increases the rate of stimulated emission, increasing y{x), Eq.|^ resulting 
in increasing separation between the relaxation oscillation side peaks and the central las¬ 
ing frequency. In all three panels of Fig.[^ the central frequency, (Uo, chosen to evaluate 
Eq. d60b is the central frequency found by the FDTD simulations. Intensity is plotted on 
a log scale in arbitrary units, rates are given in units of c/L, and the inversion and total 
number of atoms are given in SALT units of AnQ^/Tiy^^. 


general in its derivation, and can be used to predict the linewidth of semiconductor lasers given 
the appropriate form of the electric susceptibility. However, implementing an FDTD simulation 
algorithm appropriate for semiconductor gain media is challenging even in the absence of the 
effects of stimulated emission ll56l - l60l . Here, we test the N-SALT linewidth predictions using 
two-level atomic gain media; the appropriate a factor in this case was first derived by Lax 
as 121, 


Oo = 




7 ± 


(63) 


For the simulations here, we choose the atomic transition frequency almost exactly in between 
the two proximal cavity resonances, and decrease y^, thus increasing ocq. However N-SALT 
predicts a generalized a factor, a, which is sensitive to the spatial hole-burning of the gain 
medium and the non-Hermitian nature of the lasing mode, and not just to the distance of the 
lasing mode from the center of the gain curve. From ref. lIZTl in the single mode case it take the 
form: 


- _ Im[Cn] 
“ Re[Cii]’ 


(64) 


where the complex coefficients C^v determine the relaxation rate of modal fluctuations away 
from the steady-state lasing values. These coefficients are calculated from the SALT solutions 
according to: 


— 




(65) 


where (x) is the spatial profile of the /ith lasing mode, still power normalized, / {x)dx = 
1. Furthermore, it was found in Sec. 7A-B in Pick et al. Ezl, that the spatial profile of the 
first threshold lasing mode changes discontinuously as the passive cavity dielectric constant is 
increased, jumping when the hrst lasing mode switches from one passive cavity resonance to 
the next as different resonances enter and leave the bandwidth of the gain medium. Near these 
discontinuities, a large deviation between Oq and a can be observed, and we will exploit this 
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Fig. 7. Plot of the linewidth versus the output power for a two-sided dielectric slab cavity, 
n = 3.5, showing the comparison between the N-SALT linewidth prediction (green line), 
the N-SALT linewidth without an a factor (cyan line), the N-SALT linewidth using Lax’s 
a factor (blue line), and the FDTD simulation results (magenta triangles). Excellent quanti¬ 
tative agreement is seen between the FDTD simulations and the correct N-SALT linewidth 
prediction, confirming the form of the a factor derived by Pick et al. (27) . For the two-level 
gain medium used here, a„ = 18.3, Tl = 0.05, ^ = 0.01, 6 = 4 x 10^®, and Na = 10*°, 
and results in the total system having = 2.56, while ~ 0.66. Frequencies and rates 
are given in units of c/L, while the atomic values are given in SALT units of ATtO^/hyj^. 


phenomenon in the simulations below while maintaining an index of refraction similar to that 
of GaAs, using n = 3.5 for the dielectric slab cavity studied here. 

In Fig.|2]we show the results of a comparison between the N-SALT linewidth predictions 
using three different a factors, = 0 (cyan line), fa 0.66 (green line), and Uq = 2.56 (blue 
line), with direct FDTD simulation (magenta triangles). We find excellent agreement between 
the correct N-SALT linewidth calculated using a and the FDTD simulations, demonstrating 
that this is the correct form of the a factor. These simulations also verify that the Langevin 
noise model used in the FDTD simulations implicitly contains the physical effects that yield 
the a factor. While the a factor for many semiconductor lasing materials is determined ex¬ 
perimentally ED, rather than analytically, these results indicate that the physical origins of 
the phenomenon are effected by the geometry of the cavity and the spatial profile of the las¬ 
ing mode. Furthermore, this suggests that using fabrication techniques to control the index of 
refraction of semiconductor based laser cavities should allow for the engineering of different 
linewidth enhancement factors. 


5.4. two mode lasing 

A final feature of N-SALT linewidth theory is that it predicts the linewidths in the multimode 
steady state, and finds that the linewidths are not independent of one another but couple through 
gain saturation. Specifically it predicts that the linewidth of any active lasing mode is affected 
by the onset of additional lasing modes at higher pump powers. This coupling phenomenon 
occurs through a change in the a factor of each active mode at each subsequent threshold. 
Above the second lasing threshold this correction is given by 
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Fig. 8. Plot showing the modal output intensity as a function of the gain medium pump 
strength Dq, for a two-sided system consisting of two coupled dielectric cavities, n = 3, 
with different lengths, L\ = .42Lo> and ^2 = •5Lo> joined together hy a region of air, n = 1, 
with length Lair = -OSLo, where Lq is the total size of the system, and as shown in the 
schematic. This cavity has up to two active lasing modes (red and orange) for the pump 
values simulated here, and quantitative agreement is seen between the SALT simulations 
(solid lines) and noisy FDTD simulations (squares). A slight offset in the interacting thresh- 

( 2 ) 

old for the second lasing mode is seen between the two simulations, with = 0.5077, 

( 2 ) 

while Dpjjjjj = 0.5282. The inset plot shows the FDTD simulated intensity of the second 
lasing mode through its lasing threshold, first showing amplified spontaneous emission, 
then super-linear behavior at threshold, and finally linear behavior above threshold, as ex¬ 
pected. The gain medium was chosen to have cOa = 15, Y± = 0.4, ')^| = 0.01, 6 = 10^®, and 
Na = 10*®. Frequencies and rates are given in units of c/Lq, while the field quantities and 
inversion values are given in SALT units of 46^/Ti^ypy]]^ and 47t0^/hyp, respectively. 


in which single-mode N-SALT linewidth prediction from Eq. |2l and the su¬ 

perscripts R and I denote the real and imaginary components of the amplitude relaxation rates 
Cij respectively. Near threshold this analytic expression diverges and is not valid, but N-SALT 
dynamical equations can be integrated numerically to calculate the increase in the linewidth 
due to the second lasing mode M- 

To study this effect, we used two coupled dielectric cavities as shown in the schematic of 
Fig. |8l with the total system open on both ends. By using coupled cavities we create doublets 
of resonances, and by then placing the gain frequency center near one doublet we restrict our¬ 
selves to two mode lasing, which is convenient for the FDTD simulations in particular. The 
semiclassical prediction for the modal intensities as a function of the pump strength for this 
cavity calculated using SALT (solid lines) is shown in Fig.|^ and compared against the FDTD 
simulations (squares), demonstrating quantitative agreement. The inset plot shows the super- 
linear behavior through the lasing threshold observed in the FDTD simulations, as the amplified 
spontaneous emission yields a coherent lasing signal. We observe excellent quantitative agree¬ 
ment between the N-SALT prediction and the FDTD simulations for the linewidth of the first 
lasing mode, as shown in Fig.|9] On this scale, the single-mode N-SALT prediction (green) is 
very similar to the multi-mode prediction, Eq. |^(red). However, the inset of Fig. |9]shows the 
same set of comparisons through the turn on of the second lasing mode. Unfortunately, while 
there is clearly enhanced noise near the second modal threshold, there is not enough resolution 
in the FDTD data to compare accurately the single-mode and two-mode N-SALT predictions. 
There are two difficulties with the numerical comparison. First, the design of the system, and 
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Fig. 9. Comparison of the single mode N-SALT linewidth prediction (green line), two¬ 
mode N-SALT linewidth prediction (red line), and FDTD simulations (magenta triangles) 
for the first lasing mode in coupled cavity system from Fig. Inset shows a zoom in of 
the same quantities close to the interacting threshold of the second lasing mode. The two 

( 2 ) 

slightly different second mode thresholds are marked in the inset, (dashed blue line), 

12 ) 

and (dashed cyan line). While the data is too noisy, and the difference between the 

single mode and two-mode predictions too small, for the resolution of their differences, we 
do observe increased linewidth and variance in our simulations close to the threshold of the 
second lasing mode, as expected. 


the heirarchy of parameter scales that must be achieved above the floor of the spectral resolu¬ 
tion of the simulation results in a noisy signal. Second, due to discretization errors, the SALT 
and FDTD simulations give slightly different predictions for the location of the second lasing 
threshold. Thus, when plotted against the pump strength, we expect the linewidth increase in 
the FDTD simulations to occur at a slightly shifted location relative to the N-SALT results, but 
the noise makes it difficult to extract this shift. 

Although the FDTD simulations are not sensitive enough to observe the small corrections in 
the first mode linewidth due to onset of the second mode, these simulations are able to validate 
the N-SALT linewidth prediction for the second lasing mode. Figure [TO] shows a comparison 
between the N-SALT prediction and FDTD simulations for the linewidth of the second lasing 
mode as a function of the input pump strength. As noted, the offset observed between the 
two linewidths is due to the slightly different locations of the second mode threshold and if 
this difference is subtracted, as is seen in the inset of Fig. (TO] we see excellent quantitative 
agreement between the two sets of simulations. 

6. Summary 

In this work we have performed a first principles test of the N-SALT linewidth results derived 
by Pick et al. El. To do this, we used the FDTD algorithm to simulate the Maxwell-Bloch 
equations coupled to a set of Langevin noise equations, thus including the effects of spon¬ 
taneous emission. We found excellent quantitative agreement between the N-SALT linewidth 
predictions and the FDTD simulations, while finding substantial deviations from the ‘fully cor¬ 
rected’ Schawlow-Townes theory, demonstrating that the intertwining of the cavity decay rate, 
Petermann factor, incomplete inversion factor, bad-cavity correction and Henry a factor in the 
N-SALT linewidth formula is necessary and correct. This comparison was first done in a pa¬ 
rameter range in which the relaxation oscillations were weak (near the Class A boundary). 
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Fig. 10. Comparison of the single mode N-SALT linewidth prediction (green line), two¬ 
mode N-SALT linewidth prediction (red line), and FDTD simulations (magenta triangles) 
for the second lasing mode in coupled cavity system from Fig. [8] Inset shows the same 
data except with the FDTD simulations plotted at shifted pump values (cyan triangles) to 
account for the slightly different second lasing mode thresholds seen in Fig. [8] Quantitative 
agreement between the FDTD and N-SALT linewidth predictions is seen in both versions 
of the plot, but the inset demonstrates that most of the discrepancy seen in the outer plot 
is due to differences in the output power of the cavity due to the SALT simulations being 
further above threshold than the FDTD simulations for the same value of the pump Dq. 


Through comparison with the Chong and Stone linewidth theory ||25]| . we demonstrated that 
for the small, 20Xa ~ L, cavities studied here, much but not all of the improved agreement 
found by N-SALT is due to the proper treatment of the incomplete inversion factor. Next, we 
successfully demonstrated that N-SALT gives the correct output intensity spectrum including 
relaxation oscillations for Class B lasers, and correctly reproduces the side-peaks due to re¬ 
laxation oscillations. This set of simulations also verified that the side-peaks shift away from 
the center of the spectrum as the pump on the gain medium is increased. We then studied the 
different predictions for the linewidth enhancement due to the coupling between intensity and 
phase fluctuations, the a factor, and demonstrated that the N-SALT form of the a factor yields 
quantitative agreement with the FDTD simulations, while previous forms of the a factor are 
shown to disagree. This set of simulations is particularly remarkable, because in the absence of 
the N-SALT prediction for a, one might conclude that the FDTD simulations do not correctly 
capture the effects of the a factor. Instead, it is clear that the FDTD algorithm used does con¬ 
tain all of the relevant physics, and that there can be a significant difference between the various 
forms of the a factor. Finally, we demonstrated that the N-SALT theory correctly predicts the 
linewidth for multiple active lasing modes. 
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